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ABSTRACT 

Two families of models are currently considered to describe an accretion flow 
onto black holes and production of the observed X-ray radiation: (1) a standard 
cold accretion disk with a hot corona above it and (2) an outer truncated accretion 
disk with a hot semispherical inner flow. We compute spectra in the scenario with 
a hot inner flow surrounded by a truncated accretion disk covered by a hot corona 
and test the results on a sample of high redshift (z > 4) quasars observed with 
Chandra. We find that in order to reproduce the ratio of optical to X-ray fluxes 
(the a ox parameter), the optical depth of the Comptonizing plasma has to be 
rather low (r = 0.02 — 0.15 in the corona above the disk, and r = 0.15 — 0.70 
in the hot inner flow). This, together with the observed X-ray photon indices, 
implies either a high temperature in a thermal plasma (kT e = 90 — 500 keV), 
or a nonthermal electron distribution in the plasma. We put an upper limit 
on the disk truncation radius, r tr < 40i?s- The modeled accretion rate is high, 
fa > 0.2 MEdd; which may suggest that high- z radio quiet quasars are analogs of 
X-ray binaries in their high or very high state. 



Subject headings: accretion, accretion disks — galaxies: high-redshift — quasars: 
general — X-Rays 
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1. Introduction 

Modern X-ray observatories can detect the most distant objects in the universe and 
provide high quality data for quasars at redshifts as high as ~ 6 (Becker et at. 2001; Brandt et 
al. 2002b; Bechtold et al. 2003, hereafter B03), i.e. for objects born in the very early universe. 
By comparing high- and low-redshift sources, we may be able to determine whether these 
objects evolve and if there is any particular pattern that they follow during the evolution. 
This provides bases for studying the formation and evolution of the universe. 

We focus our studies on the radio quiet quasars (RQQs), which comprise ~ 85% — 90% 
of the quasar population (e.g., Stern et al. 2000). The fact that RQQs are more common, and 
thus more typical, among the quasars makes it important to understand the origin of their 
broad-band spectra. In these systems a contribution from possible additional components 
related to the radio emission should be negligible. On the other hand, RQQs are X-ray 
quieter than radio loud quasars (RLQs) in terms of higher a ox parameter — the spectral 
index of a reference power law connecting the optical and X-ray rest frame fluxes, for which 
we use the following definition: a ox = — (log/ x — log/ Q ) / (logz/ x — logz/ Q ), where f Q and / x 
are rest frame optical and X-ray fluxes at A = cjv = 2500 A and E = hu x = 2 keV (e.g., 
Zamorani et al. 1981; Worrall et al. 1987; Laor et al. 1997) — so the X-ray data of RQQs 
are of worse quality than those of RLQs. 

Understanding observed spectral features in terms of theoretical model parameters can 
provide important insights into physical processes taking place in quasars, especially about 
the origins of the hard X-ray radiation. The "standard" model so far successful in explain- 
ing optical/UV/X-ray spectra of active galactic nuclei (AGNs) (and those of X-ray binaries, 
hereafter XRBs) contains as a crucial element an optically thick accretion disk around a 
central black hole. Gravitational energy is dissipated and thermalized in the accretion disk 
(Shakura & Sunyaev 1973). It is emitted in the form of thermal radiation, giving a char- 
acteristic bump in the optical/UV band ("big blue bump") in AGN spectra (Shields 1978, 
Czerny & Elvis 1987, Siemiginowska et al. 1995, Koratkar & Blaes 1999) or soft X-rays in 
the case of XRBs (see, e.g., a review by McClintock & Remillard 2003). In the case of 
high-z RQQs, the central black hole is expected to be supermassive, with a mass of up to 
~ 1O 1O M (Rees 1984). Such a high mass is suggested by the high luminosities of high-z 
RQQs (of the order of ~ 10 48 ergs s _1 ), which should not exceed the Eddington luminosity 
of L Edd « 1.3 x 10 46 (M/1O 8 M ) erg s" 1 . 

Some of the disk photons gain energy in the inverse Compton process by scattering off 
energetic electrons in a hot, optically thin plasma. This mechanism is generally accepted to 
explain hard power-law-like spectral X-ray tails extending beyond 100 keV (see Mushotzky, 
Done & Pounds 1993 for a review; Done 2001). The spectra commonly show cut-offs beyond 
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~ 100 keV in both Seyfert galaxies and XRB (in their hard/low states), suggesting a thermal 
energy distribution of Comptonizing electrons (Gierlihski et al. 1997, Zdziarski et al. 1998). 
Observed X-ray spectra of XRBs (soft/high states) may extend to energies beyond 500 keV 
(e.g., Gierlihski et al. 1999) with no cut-off, implying a non-thermal energy distribution 
or, possibly, a hybrid (i.e., a mixture of thermal and non-thermal; Coppi 1999) electron 
distribution. 

Some of the Comptonized X-ray photons are likely to illuminate the accretion disk, 
where they can be reprocessed and mostly thermalized, thus increasing the disk flux. Since 
the disk flux provides the seed photons for the inverse Compton processes, a feedback loop 
between the cold disk and the hot plasma is created (Haardt & Maraschi 1991, 1993). 

Some of the illuminating X-ray photons will be reflected by Compton scattering off 
cold electrons and emerge back from the disk, with a characteristic spectrum (Lightman & 
White 1988; George & Fabian 1991). The spectrum is composed of (1) a broad continuum 
("reflection hump") peaking at 20-30 keV, whose shape is determined by a competition 
between photo- absorption in the soft X-ray band (E < 10 keV) and Klein-Nishina effects 
on the Compton scattering cross-section in the hard X-ray band, (2) the Fe K a fluores- 
cence/recombination line at 6.4-6.97 keV (depending on the Fe ionization state), and (3) 
the Fe-K Q edge at ~ 7.1-10 keV. The fluorescence features can be further influenced by 
relativistic effects and/or the ionization state of the reflecting medium (Fabian et al. 1989; 
Zycki & Czerny 1994, Nayakshin 2000). 

An unresolved issue is that of the geometry of the accretion flow, i.e., a geometrical 
configuration of the two phases of accreting plasma (cold, optically thick; and hot, optically 
thin). There are two general scenarios considered in the modeling: (1) the plane-parallel 
geometry, in which a hot corona is located atop the accretion disk, which is assumed to 
extend down to the innermost stable orbit (r in = 3 i?s for a Schwarzschild black hole, where 
Schwarzschild radius Rs = 2GM/c 2 ), and (2) the spherical geometry, which assumes that 
the cold disk is truncated at a certain radius r tr > rj n and surrounds a hot inner plasma flow 
within r < r tr (Shapiro et al. 1976; Poutanen et al. 1997; Esin, McClintock & Narayan 1997; 
Rozahska & Czerny 2000). In the first geometry, the hot corona sandwiching the disk may 
be continuous (e.g., Haardt & Maraschi 1991) or "patchy" i.e., composed of a number of 
active regions driven by, e.g., magnetic activity (Galeev, Rosner & Vaiana 1979). A variant 
of the spherical geometry (case 2 above) is the spherical cloud model of Collin-Souffrin et al. 
(1996), in which the cold plasma forms small clouds accreting roughly spherically onto the 
central X-ray source. 

In principle, one can distinguish between the two basic geometries by studying the X-ray 
spectrum in detail (see Done 2001 for a review). One argument is based on the amplitude 
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of the reflection component. The amplitude is defined as 

R = Q/2n, (1) 

where Q is a solid angle subtended by the reflector as seen from the source of X-rays. For 
plane-parallel geometry, R should yield values close to unity (the disk covers ~ 2ir from the 
point of view of the X-ray source), whereas for spherical geometry R < 1 would be expected. 
Furthermore, the amount of broadening of the reflection features gives an estimate of the 
inner radius of the accretion disk (e.g., Done et al. 2000). For example, the Fe Ka line would 
be narrower if the disk were truncated far away from r- m . 

Another argument in determining the geometry may be based on the X-ray photon 
index, T. A continuous, plane-parallel corona produces rather soft X-ray spectrum, T > 2, 
even when the entire gravitational energy is dissipated in the corona (Haardt & Maraschi 
1991, 1993). This is because the underlying cold disk provides a lower limit to the cooling 
flux. The spectra in the spherical geometry are harder, T < 2, because the disk flux entering 
the central sphere is strongly reduced, compared to that in the plane-parallel geometry (e.g., 
Poutanen et al. 1997). 

However, these diagnostics do not determine the geometry in a unique way, since it 
is possible to modify the continuous plane-parallel corona scenario in such a way that the 
predicted spectra are hard, with T < 2, and the amplitude of reflection is reduced, R < 1. 
This may be realized in a number of scenarios: (1) the hot plasma may form a patchy corona 
(Stern et al. 1995); (2) the cold disk may be covered by a hot ionized skin, decreasing the 
effectiveness of the reprocessing/thermalization (Nayakshin, Kazanas & Kallman 2000); or 
(3) the hot plasma may be outflowing at relativistic speeds, so that relativistic beaming makes 
the emission pattern significantly anisotropic (Beloborodov 1999a). In particular, both the 
spherical geometry scenario and models 2 and 3 above, reproduce the correlation observed 
between the X-ray spectral index and amplitude of reflection of both Seyfert galaxies and 
XRB (Zdziarski et al. 1999). 

An additional model diagnostic, which helps to break the degeneracy, is the ratio of ob- 
served optical and X-ray fluxes, the X-ray loudness, a ox . This parameter is of more use for 
AGNs than for XRBs, since the disk emission can be more easily observed in AGNs than in 
XRBs. This is especially true for high- z quasars, where the big blue bump may be redshifted 
into the observable optical band. On the other hand, the quality of X-ray data from these 
objects is usually not sufficient to accurately determine the amplitude of the reprocessed 
component (see Siebert et al. 1996, Lawson & Turner 1997, Cappi et al. 1997, Reeves & 
Turner 2000, Mineo et al. 2000, and Page et al. 2003 for attempts to fit reflected component 
and/or iron K a line to the data of quasars). 
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The "complex" models described above were applied to Seyfert galaxies and low redshift 
RQQs with high S/N (signal-to-noise ratio) data (recently, e.g., by Chiang 2002; Chiang & 
Blaes 2003; Janiuk, Czerny & Madejski 2001). So far the quasars' spectral energy distribu- 
tion, especially that of high redshift quasars, has been modeled by the thermal emission from 
an accretion disk (Sun & Malkan 1989, Bechtold et al. 1994a, 1994b, Tripp et al. 1994), with 
the addition of a Comptonizing medium explaining the X-rays (Band & Malkan 1988, Fiore 
et al. 1995). There have been no systematic studies of the parameter space of these models 
appropriate for high— z quasars, while the number of the observations has been increasing. 

In this paper we consider a geometry with a hot inner flow surrounded by a truncated 
accretion disk covered with a hot corona. We apply the model to the high redshift RQQs. 
Similar geometry was applied by Chiang (2002) and Chiang & Blaes (2003) to Seyfert galax- 
ies. Their Compton scattering medium was located within the central sphere, which was 
surrounded by an accretion disk. They extended the model of Zdziarski et al. (1999) by 
introducing viscous dissipation in the cold disk and modeled X-ray data simultaneously with 
the optical/UV data. In our model we include two Comptonizing media, the sphere and the 
corona. We explore the model parameter space and identify the model parameters suitable 
to generate the observed spectra of RQQs. The observational constraints include three quan- 
tities: the X-ray photon index, T, the X-ray loudness, a ox , and the optical/UV luminosity, 
Zuv = log iyL v ) at 2500 A in the rest frame. Optical/UV luminosities of high-z RQQs are 
of the order of Zuv ^ 46.0-47.0 (B03 and references therein). B03 and Vignali et al. (2003a, 
hereafter V03; 2003b) derived the mean ~ 2-30 keV rest frame X-ray photon index to be in 
the range ~ 1.5-2. The X-ray loudness was found to be ~ 1.7-1.8 (B03; V03; Vignali et al. 
2003b). 

The structure of the paper is as follows: In Section 2, we describe the geometry of 
accretion flow and describe the model. In Section 3, we present the results of computations, 
and discuss the constraints imposed on the parameter space by the observations. In Section 
4, we apply the spherical geometry accretion model to the high- z RQQ data of B03. Finally, 
Sections 5 and 6 contain discussion and concluding remarks, respectively. 

In a subsequent paper (Sobolewska et al. 2004, hereafter Paper II), we study the ac- 
cretion disk model covered by a non-uniform corona composed from hot clouds (a patchy 
corona) . 



- 6- 



2. Model 

We consider an accretion flow consisting of (1) a cold, optically thick accretion disk 
with (2) a hot corona above it. At a radius r tr greater than the radius of the innermost 
stable orbit, r in , the disk evaporates completely to (3) a hot inner flow. The geometry 
is schematically presented in Figure 1. The computed spectra are a superposition of the 
thermal disk radiation that escapes the hot plasma (both the corona and the inner flow) 
without being scattered, and the Comptonized component. The spectra are computed in a 
face-on view. 

Throughout the paper we assume the Schwarzschild geometry, the innermost stable 
circular orbit is located at r in = 3 Rs, where Rs = and the accretion efficiency is 

e—1/12 (unless stated otherwise). We define the Eddington accretion rate as eM-Ed&c 2 = 
1.3 x 10 46 10 £* M ergs s _1 and use the unitless quantity fn = M / M Edd to describe the accre- 
tion rate. 



2.1. The disk component 

At each radius r > r tr , we assume that the hot corona above the accretion disk dissipates 
a fraction / of the locally released gravitational energy and that the remaining energy is 
dissipated in the disk. 

The total dissipative luminosity available between the innermost stable circular orbit, 
r in , and the transition radius, r tr , is 

/■rtr 

^available = 4vr / F visc (r)rdr, (2) 

Jr in 

where -F v isc( r ) = 3< &fJ? (l — ^ e calculate the total luminosity dissipated in the hot 

inner flow as 

-^diss = ^available; (3) 

where < 8 < 1 is the radiation efficiency correction to the Keplerian efficiency within the 
inner flow. 

The cold disk (r > r tr ) is illuminated by both the hot inner flow and the hot corona. 
Part of this illuminating flux is reprocessed in the disk. Therefore, at each radius three 
contributions account for the total accretion disk soft flux: (1) viscous dissipation in the 
disk, (2) reprocessing of the hard coronal radiation, and (3) reprocessing of the hard radiation 
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originating in the inner flow (see Figure 1). The total soft disk flux is thus given by 

iWO = (l-f)F visc (r) + f V (l-a)F visc (r) + 

+ (1 - a)F inc (r), (4) 

where a and rj are the energy-integrated disk albedo and the Compton anisotropy parameter 
(Haardt & Maraschi 1991), respectively (we assume typical values of a = 0.2 and rj = 0.5; 
Haardt & Maraschi 1993), while (1 — a)F inc (r) is the flux from the inner flow intercepted by 
the disk annulus at r and F inc (r) = 16 3 ^T r / Q amax sin a cos ada ld(f> (Chiang 2002). For 

completeness, we give expressions for a max , max , and l(r,a,(f),r s ) in the Appendix. Deriva- 
tion of F inc was performed for a uniform, isotropically radiating inner flow. 



We assume that in the disk/corona region the accretion disk radiates locally as a black- 
body with a temperature T = (-F so f t /crsB) 1//4 , where <tsb is the Stefan-Boltzmann constant, 
and F so f t is given by equation (4). This local disk emission is Comptonized in the corona. The 
input spectrum for the Comptonization in the inner flow is determined as a superposition of 
local blackbodies originating in the disk at radii n, between r tr and r out , chosen arbitrarily 
to be 10 5 i?s- Each contributing blackbody is normalized to the luminosity intercepted by 
the inner flow from a given radius, L so f t (r{), given by 

Wn) = 4vr f g(r)^±rdr, (5) 

where g(r) is a geometrical factor (Chiang 2002; see also the Appendix). These contributions 
sum to the total disk luminosity intercepted by the hot inner flow, i.e., 

£ S oft = y^-^soft(n)- (6) 

n 



2.2. The Comptonized component 

In order to compute the Comptonized component, we use the empirical formulas of 
Beloborodov (1999b) that relate the amplification factor, A, the photon index, T, and the 
Compton parameter, 

y = 40 e (l + 40 e )r(l + T) (7) 

(where 9 e = kT e /m e c 2 ): 

Y= 7 -{A-l)-^\ (8) 

r = \y-V\ (9) 
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In general, the amplification factor is defined as follows: 

a= Lsoft + Ldiss , (10) 

Asoft 

where L diss stands for the total luminosity dissipated in the hot plasma and L so f t denotes the 
soft luminosity that enters the hot plasma from the accretion disk. 

We perform our analysis for electron temperatures in the hot plasma of kT e = 150 and 
500 keV. We compute the optical depth of the plasma, r, and the photon index, T, combining 
equations (7), (8), and (9). 

The relative strengths of the power-law-like high energy tail and the soft input spectrum 
escaping the hot plasma without any shift in the energy are determined by the probability, 
P T , that a soft photon is scattered to a higher energy in the Comptonization process. To 
calculate P T , we use the formula 

lnP T , , 

T - 1 = — 11) 

In A.E K ' 

where V is the photon index (r — 1 is the spectral index) and AE = 1 + 40 e (l + 48 e ) denotes 

the mean energy amplification in a Thomson collision. 

We derive formulae for the amplification factor, A, as follows. For the disk/corona 
system A, defined by equation (10), can be rewritten in terms of the fraction of the energy 
released in the corona, / (Haardt & Maraschi 1991): 

/ 

^corona = 1 + " Jf f T~^r- (12) 

1 - f[l - {l-a)rj\ 

To obtain the amplification factor in the hot inner flow, A inner , we substitute the L diss 
specified by equation (3) and the L soft given by equation (6) into the equation (10). 



2.3. Efficient, hot, spherical flow 

In the simplest version, we assume that the hot inner flow efficiency is comparable to 
the disk efficiency, i.e., 5 = 1 in equation (3), and that the hot plasma forms a sphere with 
a radius r tr . 

The model parameters are 

1. the mass of the black hole, M; 

2. the accretion rate, m; 
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3. the fraction of gravitational energy dissipated in the corona, /, (assumed constant with 
radius) ; 

4. the transition radius, r tr , at which the accretion disk evaporates to the spherical inner 
flow; and 

5. the plasma temperature, kT e (same in both Compton components). 



2.4. Low efficiency, hot, semispherical flow 

The hot plasma located within the truncation radius, r tr , may form an advection domi- 
nated flow, which was found to be several times less efficient than the accretion through the 
optically thick disk (see review by Narayan, Mahadevan & Quataert 1998). Chiang (2002) 
inferred the radiative efficiency in the Seyfert galaxy NGC 7469 to be ~ 2%-3%. Thus, we 
also consider a case with 5 < 1. In our computations we assume S = 0.12, which corresponds 
to 1% accretion efficiency In addition to the parameters listed in Section 2.3, we have now 
also the accretion efficiency correction, 5. 

In both cases (efficient and low efficiency inner flow), we check the possibility that 
the inner flow is flattened, i.e., the shape of the plasma within r tr is ellipsoidal instead 
of spheroidal. Flattened, disk-like distribution of Comptonizing plasma was suggested by 
Chiang (2002) to explain spectra of the Seyfert galaxy NGC 3516. This modification adds 
the semi-minor axis of ellipsoidal flow, r s < r tr , to the list of model parameters. 

The basic quantities that define the geometry in our approach are L^ ss and L so f t . Al- 
lowing for 5 < 1 affects both L diss and L so ft, since L soft depends on L diss through F soft ; see 
equations (5) and (4). Considering a flattened shape of the inner flow affects only L so ft, i.e. 
the soft flux intercepted by the inner flow. 



3. Results 

Based on the computed model spectra, we calculate spectral characteristics, such as 
the X-ray loudness, a ox , the X-ray photon index, T, and the optical/UV luminosity, Ijjv = 
log (yL u ) at 2500 A in the rest frame. Then we explore the model parameter space to find 
theoretical parameters compatible with the observed values of a ox , T, and I in high-z RQQs. 

Figure 2 illustrates the sensitivity of the model spectrum to different model parameters. 
As expected, the bolometric luminosity depends on the accretion rate, with almost no change 
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in a ox and T (Fig. 2b). The most dramatic spectral change is due to variations in the 
truncation radius, r tr (Fig. 2c), the electron temperature (which physically depends on the 
rate of cooling by soft photons; Fig. 2d), and the radiation efficiency of the central sphere, 
6 (Fig. 2f). 

Below, we discuss in detail the results with respect to the data of high-redshift RQQs. 
In particular, we search for values of the theoretical model parameters that give V < 2.3 and 
1.5 < a ox < 1.8. Such a choice is motivated by the 3a confidence intervals of T given by V03 
(1.7 < T<2.3), the results of Vignali et al. (2003b; T = 1.861°;^), and the B03 sample, with 
a mean photon index of T = 1.50 ± 0.15 (the error represents the 90% confidence interval). 
We choose the value of a ox m i n = 1.5 for the lower limit, motivated by the value within the 
3cr confidence interval found for RQQs in the Bright Quasar Survey (see Brandt et al. 2002a 
and references therein). The upper limit of a oxmax = 1.8 was chosen on the basis of the 
samples of B03, with a ox =1.71±0.02 (the contributions to the error of a ox come only from 
errors of the X-ray photon index and the normalization of the power law fit to the data; 
in addition, the 2 keV flux in the source frame was calculated from the 1 keV flux in the 
observer frame, assuming the photon indices found from fits, not the value of 2.2 as assumed 
in B03), and V03, with a ox =1.77±0.03 (we cite the value of a ox corrected by the authors in 
a subsequent paper, Vignali et al. 2003b). 



3.1. Efficient, Hot, Spherical Flow 

The rest frame ultraviolet luminosity, Zuv = log vL v at v — cj ^ 2500 ^ > depends on all 
model parameters except the plasma temperature, kT e . For a black hole mass of the order of 
10 10 Mq, the luminosity Zuv ~ 46.0-47.0 can be obtained for any considered accretion rate if 
the corona above the truncated disk dissipates no more than about 10% of the gravitational 
potential energy, and the remaining ~ 90% is dissipated in the disk (see Figure 3a). If the 
corona dissipates most of the energy, the UV bump becomes less pronounced, and the most 
luminous objects (with Zuv ^ 47.0) cannot be explained even with Eddington accretion rates 
unless the black hole mass is higher than 10 10 M . 

Neither the X-ray slope, T, nor the X-ray loudness, a ox , depends on the accretion rate 
for a given truncation radius and fraction of energy dissipated in the corona (see Figure 2b). 
The change in the accretion rate influences only the normalization of the spectrum. 



The X-ray photon index, T, depends on the strength of the corona, /, and the disk 
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truncation radius, r tr (Fig. 2a and Fig. 2c). The strength of the corona needs to be fixed 
at ~ 10% in order to match the observed optical/UV luminosity (as stated above), so the 
only parameter that determines T is the disk truncation radius, r tr . If the truncation radius 
increases, the amplification factor of the hot inner flow, v4 inncr , increases, and a harder X-ray 
spectrum is produced. For 4i?s < r tr < 100 R$, the model gives an X-ray photon index 
ranging from ~ 2.3 to ~ 1.6. 

The X-ray loudness, a ox , depends on the strength of the corona above the disk, /, the 
truncation radius, r tr , and the plasma temperature, kT e . The strength of the corona, /, is 
constrained by the optical/UV luminosity. An increase in the disk truncation radius, r tr , 
causes the X-ray loudness to decrease. This is because an accretion disk truncated farther 
away from the black hole produces a less pronounced optical/UV bump, and Zuv decreases. 
Simultaneously, the X-ray spectrum hardens (as described above), and a higher 2 keV flux 
is produced (see Figure 2c). 

Variations of the electron temperature influence significantly the X-ray loudness, a ox , 
but do not influence the X-ray slope, T (Fig. 2d). The maximum value of the X-ray loudness 
for a temperature of 150 keV is about 1.7. For a higher temperature, kT e = 500 keV, the 
maximum X-ray loudness, a ox ^ max , could be as high as ~ 2.2. 

The lack of the dependence of the X-ray slope on the plasma temperature follows from 
the fact that in the model T is determined by the amplification factor; i.e., the X-ray slope is 
only geometry dependent. The sensitivity of the X-ray loudness to the plasma temperature 
can be explained by the fact that in thermal plasmas spectral characteristics, such as the 
X-ray photon index, the optical depth, and the electron temperature, are not independent. 
An increase in the electron temperature is compensated for by a decrease in the optical 
depth, resulting in approximately the same X-ray slope. The higher electron temperature 
implies the higher energy of the spectral cut-off. The total energy must, however, be con- 
served. Therefore, the normalization of the Comptonized component drops, resulting in the 
X-ray quieter spectra. Plasma temperature variations only weakly influence Zuv ; through 
the dependence on the amount of the soft radiation that escapes the plasma without being 
Comptonized on kT e (see Equation (11)). However, this effect is not significant enough to 
influence the results. 

Figure 4 summarizes the above results and shows two observationally interesting re- 
gions, 1.5 < a ox < 1.8 and 1.7 < T < 2.3, in the truncation radius-accretion rate plane. The 
computations were performed for two electron temperatures, kT e = 150 and 500 keV, and 
two values of the fraction of energy dissipated in the corona, / = 0.1 and / = 0.6. For an 
electron temperature of 150 keV and / = 0.1, the two regions overlap if the transition radius 
is of the order of several Rs- The spectra are relatively soft and X-ray loud, with T ~ 2.0-2.3 
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and ct ox ~ 1.5 — 1.7 (Figure 4a). For higher electron temperatures, the overlapping region 
exists at higher truncation radii (r tr ~ 7 — 15 Rs for kT e = 500 keV). In addition, harder 
(T < 2.0) and X-ray quieter (a ox ~1.7-1.8) spectra can be obtained (Figure 4b). A stronger 
corona (/ = 0.6) requires a higher electron temperature (kT e > 200 keV) to produce spectra 
with T ~ 2.0 — 2.3 and a ox ~ 1.5 — 1.7. For kT c > 400 keV, the model additionally gives 
spectra with a ox ~ 1.7 — 1.8 and T ~ 2.0 — 2.3, or a ox ~ 1.5 — 1.7 and T ~ 1.7 — 2.0 (Figure 
4c). However, in case of / = 0.6, the optical/UV luminosity is Zuv < 46.2. 



3.2. Low-Efficiency, Hot, Semispherical Flow 

We compute a ox , T, and /uv in the low efficiency hot inner flow model (Section 2.4). 
All the spectral trends with the strength of the corona, the accretion rate, the truncation 
radius, and the electron temperature are qualitatively the same as for the high efficiency 
hot flow. As in Section 3.1, a weak corona above the disk is required in order to fit the 
optical/UV luminosity. The dependence of Zuv on the truncation radius and the accretion 
rate for / = 0.1 is shown in Figure 3b. 

Figure 2f shows dramatic changes in the X-ray slope, T, and X-ray loudness, a ox , with 
reduced efficiency of the central hot flow. For a given truncation radius, r tr , less heat is 
being supplied to the electrons (some of the energy is advected); i.e., the total luminosity 
of the hot flow decreases, Ldiss(5 < 1) < Ldi SS (S =1). The soft disk luminosity intercepted 
by the hot flow, L so f t , also decreases, because it depends on the hot flow luminosity through 
F inc (see equations (5) and (4)). However, L diss decreases faster than L so{t , since the latter 
contains contribution from the viscous dissipation in the outer disk. Thus, lowering the hot 
flow efficiency causes a decrease in the amplification factor of the hot flow, A nner , and leads 
to a softer X-ray spectrum. In addition, decreasing L so f t lowers the normalization of the soft 
input spectrum for Comptonization in the hot flow, which results in a lower normalization 
of the X-ray continuum and X-ray quieter spectra. This can be clearly seen in Figure 2f for 
the plasma temperature of 500 keV. 

To summarize, for a given truncation radius, r tr , lowering the efficiency of the hot inner 
flow leads to a softer X-ray slope and X-ray quieter spectra. This result is illustrated in 
Figure 5. To produce spectra with the same T and a ox as in the case of the efficient flow 
(compare Figure 4), a sphere with a larger radius is required. As in Section 3.1, a plasma 
temperature higher than ~ 150 keV and a fraction of energy dissipated in the corona smaller 
than ~ 10% are implied by the data. However, the transition radius now falls into the range 
of ~ 10 — 70 R$, depending on the strength of the corona, the plasma temperature, and the 
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accretion rate. The hardest, T ~ 1.7 — 2.0, and X-ray quietest, a ox ~ 1.7 — 1.8, spectra can 
be obtain for higher electron temperatures. Figure 5b shows the case of kT e = 500 keV with 
the solution for r tr ~ 30 — 40 Rs- 

In the second modification to the efficient flow, we allow for a flattening of the central 
sphere, e.g., an ellipsoidal flow. The amplification factor now depends strongly not only on 
r tr , but also on the semi-minor axis of the ellipsoidal flow, r s , and on the ratio r s /r tr . For 
a given transition radius there exists the following trend: the smaller the r s , the harder the 
spectra (see Figure 2e). This is because the luminosity, L so ft, drops, resulting in an increase 
of the hot flow amplification factor, A- mner (Ldi SS remains constant). This geometry may 
produce spectra as hard as V = 1.5, but with a ox significantly smaller than that observed 
in high-z RQQs if the hot inner flow is efficient, 5 = 1 (Section 2.3). If the radiative 
efficiency is reduced, 5 < 1 (Section 2.4), allowing for ellipsoidal central flow results in a 
shift of the regions indicated in Figure 5 toward slightly lower transition radii (i.e., it is an 
intermediate case between efficient spherical inner flow and efficient flattened inner flow). 
Thus, we consider r tr ~ 40 Rs to be an upper limit of the truncation radius in the model 
that allows us to obtain the X-ray quietest spectra with a ox ~ 1.7 — 1.8. 

We note that the objects with very hard photon indices, T < 1.5, cannot be explain by 
the model, because the X-ray loudness of such spectra is always lower than a ox = 1.5-1.8. 



4. Application to the data 

We apply the model to eight high redshift RQQ from the B03 sample. We choose only 
the seven sources with observed photon index T > 1.5. We also include the highest redshift 
quasar, at z=6.28. We do not perform a formal x 2 fitting to determine the best-fit parameters 
and their confidence limits, but rather demonstrate the parameter values that approximately 
reproduce the spectra. With the available data the results are not unique, and a number 
of solutions are possible. The computed spectra are shown in Figures 6 and 7. The fit 
parameters are listed in Tables 1 (efficient inner flow) and 2 (low efficiency inner flow). The 
objects were observed by Chandra in the 0.3-6.5 keV range. The 1450 A rest frame points 
were taken from the literature, as described in detail in B03. The ratio of the optical/UV 
and X-ray fluxes is characterized by the X-ray loudness. In the case of high-z quasars, the 
2500 A flux used to calculate a ox must be extrapolated from the 1450 A flux known from 
observations, which involves a knowledge of the optical/UV spectral index (/„ ~ u avv ). To 
derive general properties of a theoretical model or an observational sample, a mean value 
suggested by observations can be used (as in Section 3 or in B03, where auv — —0.3 was 
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adopted; Kuhn et al. 2001)). However, it does not apply to modeling the spectra of particular 
objects. Thus, in this section we use the flux at 1450 A to normalize the spectrum, and we 
list values of the spectral index between rest frame 1450 and 2500 A, «uv, and the X-ray 
loudness, a ox , computed from the model spectra in Tables 1 and 2. The optical and X-ray 
observations were not simultaneous. However, Giveon et al. (1999) argue that the optical 
variability anti-correlates with the luminosity of AGNs. Given that the high-z RQQs are 
very luminous objects, one would expect optical variability at a relatively low level. 

Matching the observed optical/UV luminosity requires weak coronae above the accretion 
disk, dissipating no more than ~ 10% of the gravitational potential energy. In addition, high 
masses of the central black holes should be considered, M ~ (0.5 — 1.4) x 1O 1O M . Our 
modeling of BRIO 103+0032 illustrates the following trend: an increase of the black hole 
mass (with the X-ray slope kept constant) forces the accretion rate to decrease in order to fit 
the optical/UV luminosity. In addition, the electron temperature drops (the plasma optical 
depth increases) to match the X-ray emission. 

For the radiatively efficient spherical inner flow (Section 2.3; 5 = 1 and r s = r tr ), the 
truncation radius of the accretion disk is typically 5-11 Rs- This is constrained mainly 
by the observed X-ray slope. In general, the electron temperatures of ~ 108-500 keV are 
required in order to match the observed X-ray loudness. The X-ray quietest objects in our 
sample (a ox ~ 1.7 - 1.8; PSS1057+4555, SDS1030+0524, and SDS1204-0021) require very 
low optical depth of the plasma (r < 0.05), which brings the temperature to unrealistically 
high values of 700-1400 keV. Thus, models with 5 = 1 are not acceptable, at least for 
PSS1057+4555, SDS1030+0524, and SDS1204-0021. 

Considering a flattened hot inner flow allows for slightly lower electron temperatures. 
However, in the case of objects with a ox ~ 1.7-1.8, the electron temperature is still high, 
ranging from 410 to 800 keV. In the case of the two objects with T > 2.3, there is no need to 
introduce this modification, since the flattening of the hot flow would cause the hardening 
of the spectrum. These soft objects can be modeled with a very small hot inner flow (with 
r tr ~ 5 Rs) or even with a standard plane-parallel configuration in which the accretion disk 
covered with a hot corona extends down to the last stable orbit, r in (in Table 1, fits with 
ftr = r s = ri n ). In addition, the corona may be stronger and can dissipate up to 70% of the 
energy (see fits for BRI1033-0327 and PSS1317+3531 in Table 1). 

Assuming the low efficiency spherical flow (Section 2.4; 5 < 1, r tr = r s ) results in 
an increase of the disk truncation radius needed to fit the observed X-ray slope. Now 
rtr ~ 9 — 35i?s (see Table 2). The electron temperature needed to match the observed 
X-ray loudness is kT c ~ 90-345 keV. However, the objects with the hardest X-ray spectra 
(SDS1030+0524 and SDS1204-0021) cannot be modelled using low efficiency spherical flow 
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geometry. Reproducing their X-ray slope (within a la confidence interval) requires the disk 
to be truncated at 60 — 80 R$. With this large an r tr , the black hole mass would have to be 
larger than 2 x 10 10 M & in order to fit the optical/UV luminosity, even assuming accretion at 
the Eddington rate. Even then the fit is not acceptable, because of the optical/UV spectral 
index (we discuss this issue in more detail below). 

We also perform fits in which we allow for both low efficiency and flattening of the 
hot inner flow. In this approach, all the objects with T < 2.3 can be modeled. The disk 
truncation radius is r tr ~ 11-20 R$, the semi-minor axis of the inner flow is r s ~ 4-10 R$, 
and the electron temperature ranges from ~ 195 to ~ 500 keV. 

The model covers a wide range of the optical/UV spectral index. For the efficient inner 
flow, we get —0.6 < auv < —0.3 and for low efficiency inner flow, —1.5 < «uv < —0.3. 
These values are in agreement with observations that suggest a significant spread in ayy] 
e.g., Fan et al. (2001) derive «uv = —0.79 ±0.34 (la), and Pentericci et al. (2003) report on 
auv = —0.57 ± 0.33 (la) for high-z quasars (3.6 < z < 5). Knowing the exact value of «uv 
for individual sources would reduce the number of model parameters and make the fitting 
procedure unique. 

We stress that the results for the plasma temperature and optical depth come from the 
constraint placed by the ratio of optical/UV to X-ray emission. The values of X-ray loudness 
we obtain from fits differ by up to 16% from those given by B03. This difference comes from 
the adoption of ctuv — —0.3 and T = 2.2 by B03 when calculating a ox . 

Our modeling implies a rather high accretion rate, M ~ 0.2-1.4 MEdd- This would sug- 
gest that the high- z RQQs may correspond to high or very high states of black hole binaries 
(see the comparison of AGN and XRB luminosity states in Czerny 2003). However, the 
accretion disk in the high or very high state of an XRB is expected to extend to the last 
stable orbit (e.g., Esin, McClintock, & Narayan 1997). 



5. Discussion 

We have computed optical/UV/X-ray spectra from a two-component accretion flow in 
the geometry with a hot, semispherical inner flow enclosed by a truncated accretion disk 
covered with a hot corona. We assumed that the soft disk blackbody-like radiation was 
Comptonized in the hot plasma and that the electron energy distribution in the hot plasma 
was thermal. We investigated the values of model parameters for which the model can 
account for spectra observed in high-redshift RQQs. Such spectra can be characterized by 
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three observables: (1) the luminosity at 1450 A in the rest frame of the source, Zuv = 
log [yLy] > 46.0-47.0, (2) the X-ray loudness, a ox ~ 1.5-1.8, and (3) the X-ray photon 
index, T ~ 1.7-2.3. Finally, we applied the model to the RQQ data of B03. 

The high luminosities of objects in our sample require not only rather large black hole 
masses, ~ 1O 1O M , but also high accretion rates, above 0.2MEdd- 

We found that some objects require the radiative efficiency of the hot inner flow to 
be reduced compared to that of the standard Keplerian disk. A significant fraction of the 
energy is then advected inwards, rather than radiated (i.e., an ADAF-like flow). The hardest 
(T < 2.0) and X-ray quietest (a ox ~ 1.7-1.8) spectra require a truncation radius up to 
r tr ~ 30-40 R$ if the inner flow is spherical. If the inner flow is additionally flattened, the 
truncation radius is smaller, yielding ~ 15-30 R$. 

These values of the truncation radius are comparable to those obtained by Janiuk et al. 
(2004), who studied a time evolution of unstable accretion disks with possible disk evapora- 
tion. Janiuk et al. considered the effects of the MHD turbulence on the viscosity during the 
evolution of the thermal-viscous ionization instability in the standard a-accretion disks, in 
order to explain intermittent activity in AGNs. They considered the possibility of accretion 
through a truncated disk whose inner part is evaporated to an ADAF-like flow. They found 
that the evaporation is also important during the outbursts, even at relatively high accretion 
rates of 0.1 MEdd- The disk then truncates at ~ 25 i?g (for a black hole mass of 10 s M & and 
viscosity parameter a = 0.1). We note, however, that their value for the disk truncation 
radius depends strongly on the adopted value of the viscosity parameter. 

We were not able to model objects with very hard X-ray spectra (r = 0.33-1.44 in the 
B03 sample). Such hard spectra can be produces within the model for the accretion disk 
truncated far away from the black hole (r tr >> 10 2 i?s)- However, the optical/UV part of 
the spectrum is then suppressed and the model a ox value is much smaller than that observed. 

The truncation radius can, in principle, be tested using the X-ray reflection component. 
In the considered geometry, most of the X-ray reflection comes from disk radii close to r tr , 
so the sharp spectral features (Fe Ka line and edge) should be broadened by the relativistic 
effects (Fabian et al. 1989; Laor 1991). The amount of reflection is given by the amplitude 
of the reflected component, R. In our approach, the accretion disk may reflect hard X-rays 
emitted by both the hot corona and the hot inner flow. The amplitude of reflection is thus 
given in our model by 

R=^, (13) 

-^hard 

where L inc is the luminosity emitted by the hot plasma (hot inner flow + hot corona) toward 
the accretion disk and Lhard is the total luminosity of the hot plasma (hot inner flow + hot 
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corona) as seen by a distant observer, i.e., 




oc 




inc 



sphere, inc 



corona, inc 



rF inc (r)dr + 2n 



rF visc (r)r]fdr 



(14) 



(we take 2ir, since we consider a top view scenario in which the observer sees just one surface 
of the disk) , and 



The reflection amplitude of only the coronal radiation yields R — 1, if one assumes isotropic 
Comptonization with rj = 0.5 (Haardt & Maraschi 1991). However, in our geometry the 
dominant contribution comes from the reflection of the inner flow radiation by the disk. 
Thus, the typical values of the amplitude of reflection are R = 0.2-0.6 (the efficient, spherical 
inner flow) and R = 0.15-0.35 (low-efficiency, semispherical inner flow). The lower value of 
the reflection amplitude in the case of the low-efficiency, semispherical inner flow is caused 
by the higher value of the disk truncation radius needed to fit the data (see Table 2). The 
higher truncation radius (i.e., the accretion disk is truncated farther away from the black 
hole) means a smaller contribution from the reflection of coronal radiation. We note that 
if the R — T correlation observed in XRBs and Seyfert galaxies (Zdziarski, Lubihski, Smith 
1999) also holds for high-z RQQs, this model may account for it. In the case of objects 
with high accretion rates, the reflecting material will be highly ionized. As a result, the 
reflected component at energies smaller than the fluorescent iron line energy will be flatter, 
as compared to neutral reflection. Given the small amplitude of the reflected component 
calculated above, we estimate that the change of the X-ray slope due to reflection will be 
not bigger than 0.1-0.2, which is well within observational errors. Tests involving studies of 
the reflected continuum and iron fluorescence line require X-ray data of good quality, with 
the uncertainties not exceeding 1%, which may be available in the future observations with, 
for example, ASTRO-E2 and the Constellation- X mission. 

Results of our modeling imply that the temperature of Comptonizing plasma should be 
rather high, up to 500 keV. This constraint comes from the required low optical depth of the 
Comptonizing plasma needed to reproduce the relatively low normalization of hard X-rays 
compared to that of the disk component (i.e., the high value of a ox ). In Figure 8, we present 
the range of the amplification factor that we obtain for the Comptonization in the efficient, 
hot, spherical inner flow (A = 2.5-60, lower axis) and the corresponding range of the optical 
depths (r = 0.3-1.0; upper axis) for the electron temperature 150 keV. Similar dependencies 
are shown for the Comptonization in the corona (A = 1.0-3.5; r = 0.05-0.37). Obviously, 
the higher the electron temperature, the lower the optical depth, because the X-ray spectral 
slope is fixed by the data. 
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Fits to the spectra of accreting black hole XRBs in the hard state and local Seyfert 
galaxies give electron temperatures of 50-150 keV and optical depths of the order of r ~ 1 
(Gierlihski et al. 1997; Zdziarski et al. 1998; Zdziarski, Poutanen & Johnson 2000; Done 
et al. 2000). In particular, Zdziarski et al. (2000) analyzed the average spectrum of 17 
radio quiet Seyfert 1 galaxies in the 50-500 keV range observed with OSSE. They chose the 
spherical geometry with an input of soft photons at the temperature of 10 eV in the center 
of the sphere. They used the compps model in the Xspec package to fit the Comptonized 
continuum together with the reflection feature. Thus, the high energy cutoff in the model 
spectrum, which determines the plasma temperature as a physical parameter, was directly 
constrained by the data points. Their derived plasma temperatures for Seyfert galaxies were 
much lower than the temperatures required for high redshift quasars in our modelling. 

On the other hand, Chiang (2002) and Chiang & Blaes (2003) performed analysis of 
three Seyfert galaxies that was very similar to the analysis in the present paper. By fitting 
the spectra including both the optical/UV and the 2-20 keV continuum, they also obtained 
rather high plasma temperatures (sometimes above 500 keV) and low optical depths. The 
difference between their approach and the approach of Zdziarski, Poutanen & Johnson (2000) 
was modeling the X-ray data together with optical/UV data. In addition, Chiang (2002) 
and Chiang & Blaes (2003) attempted to additionally constrain their model by inferring 
black hole masses and disk truncation radii from other observations (reverberation mapping, 
width of the Fe Ka line, etc.). They did not use data above 30 keV, so the hard X-ray cutoff 
was a free parameter in their model, as it was in ours. 

To qualitatively compare high- and low-z objects, we perform exemplary fits to the 
Seyfert galaxy NGC 7469 (analyzed by Chiang 2002) and the low-z (z=0.206) quasar PG 
0947+396 from a sample described in Laor et al. (1997). For NGC 7469 we obtain a good 
fit for kT e = 110 keV and r s = 0.8 (in agreement with findings of Zdziarski et al. 2000), a 
truncation radius of r tr = 27_Rg, an accretion rate of 0.17 times the Eddington rate, and a 
mass of 5xl0 7 M . The inner flow efficiency is low, comparable to that of Chiang (2002). In 
the case of PG 0947+396, we require the accretion disk to extend to the last stable orbit, the 
corona to dissipate 40% of the energy, the plasma temperature kT e to be 230 keV, and the 
optical depth r c to be 0.1. The mass estimate is 5xlO 8 M , and the accretion rate accounts 
for 0.7 times the Eddington rate. These results suggest that low-z objects can be modeled 
with lower masses and different plasma parameters (lower temperatures and higher optical 
depths). However, a satisfactory comparison between low- and high- z samples cannot be 
made at this stage. First, the X-ray data do not correspond to the same rest frame energies: 
2-30 keV for high-z quasars, 2-10 keV for NGC 7469 (Nandra 2000), and 0.3-2 keV for PG 
0947+396 (Laor et al. 1997). Second, the low-z objects require very detailed data analysis 
in the 0.3-10 keV X-ray range to model the soft X-ray excess, and above 10 keV to model 
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hard X-rays. In principle, our model is able to explain the soft excess in terms of two 
Comptonizations, since it includes two Comptonizing media: a corona above the accretion 
disk and an inner flow, unlike, e.g., the models of Chiang (2002) and Chiang & Blaes (2003). 
However, in the recent paper of Gierlihski & Done (2004), the authors model the 0.3-20 keV 
range with only one Comptonization component and a complex absorption feature and point 
out that the soft X-ray excess may be an artifact of such an absorption. The present data 
of high- z quasars do not allow a determination of whether their rest-frame spectra show the 
soft X-ray excess or not, because the component is redshifted below the observable band. 
On the other hand, some absorption features may become important in the future, higher 
quality X-ray data. 

Modeling broad-band quasar spectra including both accretion disk emission and the 
Comptonized hard X-rays provides a difficult theoretical challenge. Reproducing the rela- 
tive normalization of the two main emission components requires rather detailed knowledge 
of the geometry and emission mechanisms and patterns. Moreover, the high accretion rate 
inferred from our modeling implies that the studied objects may correspond to high or very 
high states of galactic black hole binaries. In addition, the resulting high electron temper- 
atures and low optical depths in the scenario with a thermal plasma may suggest that the 
Comptonization actually takes place in non-thermal or hybrid plasmas, as is seen in XRBs 
in their high state (i.e., Gierlihski et al. 1999). Many detailed studies of XRBs reveal com- 
plex, multi-component spectra (see, e.g., Zycki, Done & Smith 2001; Wilson & Done 2001; 
Gierlihski & Done 2003). The spectra contain not only simple disk emission with a hard 
Comptonized component, but also additional Comptonization of the disk emission, mod- 
ification due to X-ray irradiation of the external regions of the disk, and the reprocessed 
component, usually ionized and relativistically smeared. These additional components are 
also detected in AGN spectra (e.g., PG 1211+143; Janiuk, Czerny & Madejski 2001), and 
will surely have to be incorporated into future models (e.g., Czerny et al. 2003). 



6. Conclusions 

In this paper we explored a parameter space of the semi-spherical accretion models 
applicable to high- z RQQs. Our results show the following: 

1. The observed optical/UV luminosity suggests that high redshift quasars host supermassive 
black holes with masses of the order of 10 10 M Q . The modeling results in accretion rates higher 
than ~ 0.2 times the Eddington rate. 
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2. The high accretion rates may suggest similarity between high-z RQQs and XRBs in the 
high or very high state. If high- z RQQs correspond to high or very high states of X-ray 
binaries, the high values of plasma temperature (up to 500 keV) and low optical depths 
might suggest instead that Comptonization takes place in a non-thermal or hybrid plasma. 

3. The observed X-ray loudness puts strong constraints on the modeled plasma temperature 
and optical depth. Chandra or XMM-Newton observations with better S /N are needed to give 
tighter constraints on the X-ray slope. Future mission (such as Constellation- X) can provide 
information about the high-energy cut-off and the reflected component that is necessary to 
uniquely determine geometry of the accretion flow. 

4. Precise observations constraining the optical/UV spectral index for individual sources 
would help to reduce the number of model parameters and perform more accurate fits. 
However, the intrinsic reddening of the quasar spectrum is the main observational uncertainty 
in this spectral region. 

5. With the present data, a number of theoretical descriptions are possible. In Paper II, we 
focus on another interesting possibility, in which the accretion disk extends to the last stable 
orbit and the hot Comptonizing corona is formed by clouds above the disk (the "patchy" 
corona geometry). 

We thank Martin Elvis, Paul Green, Jill Bechtold, and Jim Chiang for many discussions. 
We also thank the anonymous referee for detailed reading of the manuscript and helpful 
comments. This research is funded in part by NASA contract NAS8-39073. Partial support 
for this work was provided by the National Aeronautics and Space Administration through 
Chandra Awards Number G01-2117B and G02-3148A, issued by the Chandra X-Ray Center, 
which is operated by the Smithsonian Astrophysical Observatory for and on behalf of NASA, 
under contract NAS8-39073. M.A.S. and P.T.Z. were partially supported by the grant PBZ- 
KBN-054/P03/2001 and KBN project number 2P03D00322. M.A.S. acknowledges support 
from the Smithsonian Institution Pre-doctoral Fellowship program. 



For completeness, we give here expressions for g(r), a max , max , and Z(r, a, <f>, r s ) (see 
Chiang 2002): 



A. Appendix 
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Fig. 1. — Hot inner flow geometry. The accretion disk evaporates at a radius r tr > r- m and 
forms a hot, semispherical inner flow. The total luminosity of the hot inner flow is L^ iss . 
The truncated disk is covered with a hot corona at radii r > r tr . At each radius r > r tr , 
three components contribute to the soft disk flux -F so f t (r): a viscous dissipation in the disk; 
reprocessing of the hard coronal radiation; and reprocessing of the hard radiation from the 
inner flow, F inc (r). The spectrum is composed of (1) the soft disk blackbody-like radiation 
that escapes plasma without being scattered, (2) the component Comptonized in the hot 
corona, and (3) the component Comptonized in the hot inner flow. 
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Fig. 2. — Spectra computed in geometry with the hot inner flow. Unless stated differently, 
the fixed parameters account for M=10 10 M , rh— 0.1, /=0.5, r tr =7 Rs, r s = r tr , /cT c = 150 
keV, and 5=1. The panels show the dependence of the following model parameters: (a) 
the fraction of energy dissipated in the clouds, f=0.1 (solid curve), 0.5 (dashed curve), and 
0.9 (dotted curve); (b) the accretion rate, m=0.1 (solid curve), 0.5 (dashed curve), and 0.9 
(dotted curve); (c) the transition radius, r tr = 5i?s (solid curve), 10Rs (dashed curve), and 
20i?s (dotted curve) Rs; (d) the electron temperature, kT e =150 (solid curve), 350 (dashed 
curve), and 500 (dotted curve) keV; (e) the semiminor axis of the ellipsoidal flow, r s = 25i?s 
(solid curve), 10i?s (dashed curve), 5Rs (dotted curve), with r tr = 25i?s; (f) the hot inner flow 
efficiency correction, 5 = 1 (solid curves), and 5 < 1 (dotted curves), for electron temperature 
of 150 (upper curves) and 500 keV (lower curves). 
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Fig. 3. — Dependence of the optical/UV luminosity on the disk truncation radius, r tr , and 
the accretion rate, rh, for (a) 6 = 1 and (b) 5 < 1 (low-efficiency flow). The shaded regions 
produce spectra with UV luminosity 46. 0< Zuv <47.0. The solid curves enclose the part of 
parameter space with Zuv > 47.0. The computations were done for a black hole mass of 10 10 
M Q , the fraction of energy released in the corona of / = 0.1, and the electron temperature 
of kT e =150 keV. The result does not depend on the electron temperature. If the corona 
becomes stronger (/ increases), the region with Z uv > 47.0 disappears from the plot, and 
the most luminous objects cannot be fitted unless M>10 10 M Q . 
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Fig. 4. — Results for geometry with efficient, hot, spherical inner flow (5 = 1, r s = r tr ; 
Section 2.3), showing the overlapping of regions yielding 1.5 < a ox < 1.8 and 1.7 < T < 2.3. 
The solid curve encloses the part of the parameter space with the hardest spectra (r < 2). 
The dashed curve encloses the region with the X-ray quietest spectra (a ox > 1.7). The 
contours are calculated for kT e = 150 keV (a) and 500 keV (b and c), and a fraction of 
energy dissipated in the corona of / = 0.1 (a and b), and / = 0.6 (c). The overlapping exists 
for r tr < 15i?s- For electron temperatures < 150 keV, the model can produce relatively soft 
and X-ray loud spectra with T ~ 2.0 — 2.3, and a ox ~ 1.5 — 1.7 (a). For higher electron 
temperatures, harder (r < 2.0) and X-ray quieter (a ox ~ 1.7 — 1.8) spectra can be obtained 
(b). For a stronger corona and high electron temperature, the model gives either hard or 
X-ray quiet spectra (c). The computations were done for a black hole mass of 1O 1O M . 
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Table 1. Modeling the SED of high- z RQQs. Efficient Inner Flow. 
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b Accretion rate in Eddington units, Mgdd 
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Fig. 5. — Results for geometry with low-efficiency, hot, spherical inner flow (5 < 1, r s = r tr ; 
Section 2.4), showing the overlapping of regions with 1.5 < a ox < 1.8 and 1.7 < T < 2.3. 
The overall trends are similar to those illustrated in Figure 4. However, the truncation 
radius is now ~ 10 — 70 Rs, depending on the other parameters. The hardest and X-ray 
quietest spectra are produced for r tr ~ 30 — 40i?s and high electron temperatures (b). The 
computations were done for a black hole mass of 10 10 M . 
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Fig. 6. — Modeling the objects from the B03 sample. Left Hot semispherical flow (5=1). 
Right: Low-efficiency semispherical flow (5 < 1). The curves show the following cases: 
?"tr = r s {solid, long-dashed, dot-dashed curves), r tr > r s (dashed and dotted curves), f > 0.1 
(dot- dashed curves), and plane-parallel geometry (dash-long- dashed curves). Fit parameters 
are listed in Tables 1 and 2. 
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Fig. 7. — Modeling the objects from the B03 sample. Left: Hot semispherical flow (6=1). 
Right: Low-efficiency semispherical flow (S < 1). The curves show the following cases: r tr =r s 
(solid, long-dashed, and dot-long- dashed curves), r tr > r s (dashed and dotted curves), f > 0.1 
(dot-dashed curves), and plane-parallel geometry (dash-long- dashed curves). Fit parameters 
are listed in Tables 1 and 2. 
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Fig. 8. — Hot plasma parameters in the model with no reduction in efficiency (5 = 1), 
and spherical inner flow (r tr = r s ) for (a) r < r tr (the inner flow) and (b) r > r tr (the 
corona). The thick solid curves show the dependence of the Compton parameter, y, on the 
amplification factor, A, obtained by combining eqs. (8) and (9). Vertical dotted lines (at [a] 
A = 2.5 and A = 60 and [b] A = 1 and A = 3.5) show the range of A possible to obtain 
from the model. The horizontal dashed lines show the corresponding range of y. Solid lines 
illustrate the relation between the Compton parameter, y, the optical depth, r (upper axis), 
and the electron temperature, kT e (eq. [7]). The range of y marked by the horizontal lines, 
corresponds to the range of r for a given kT e . For illustration, the range of r for kT e = 150 
keV is shown by dash-dotted lines. 
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c Fraction of gravitational energy dissipated in the corona above the disk 

d Disk truncation radius, r tr , and the semi-minor axis of the ellipsoidal flow, r s , in Schwarzschild units, _Rg 
c Plasma temperature in keV 

'Optical depth of corona, r c , and semispherical flow, T a 

s X-ray photon index allowed to vary within the lcr confidence interval (we do not list £tjv from fits since 
its value was fixed at that observed) 

h UV index between 1450 and 2500 A in the rest frame 

'X-ray loudness 
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Table 2. Modeling the SED of high-z RQQs: Low-Efficiency Inner Flow 
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a Black hole mass in Solar masses, Mq 
b Accretion rate in Eddington units, Aij;dd 

c Fraction of gravitational energy dissipated in the corona above the disk 

d Disk truncation radius, rtr, and the semi-minor axis of the ellipsoidal flow, r B , in Schwarzschild units, ijg 
c Plasma temperature in keV 

'Optical depth of corona, t c , and semi-spherical flow, t s 

s X-ray photon index allowed to vary within the lcr confidence interval (we do not list £tjv fr° m fits since 
its value was fixed at that observed) 

h UV index between 1450 and 2500 A in the rest frame 

'X-ray loudness 



